Draft version March 2, 2013 

Preprint typeset using I^T^X style cmulateapj v. 5/2/11 
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ABSTRACT 

Using the Millennium-II Simulation dark matter sub-halo merger histories, we created mock catalogs 
of Lyman Alpha Emitting (LAE) galaxies at z = 3.1 to study the properties of their descendants. 
Several models were created by selecting the sub-halos to match the number density and typical dark 
matter mass determined from observations of these galaxies. We used mass-based and age-based 
selection criteria to study their effects on descendant populations at z ~ 2, 1 and 0. For the models 
that best represent LAEs at z = 3.1, the z = descendants have a median dark matter halo mass of 
10 12 - 7 Mq, with a wide scatter in masses (50% between 10 118 and 10 13 7 M ). Our study differentiated 
between central and satellite sub-halos and found that ~ 55% of z = descendants are central sub- 
halos with MModian ~ 10 12 . This confirms that central z = descendants of z = 3.1 LAEs have halo 
masses typical of L* type galaxies. The satellite sub-halos reside in group/cluster environments with 
dark matter masses around 10 Mq. The median descendant mass is robust to various methods of 
age determination, but it could vary by a factor of 5 due to current observational uncertainties in the 
clustering of LAEs used to determine their typical z — 3.1 dark matter mass. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: high-redshift - large-scale 
structure of universe 



1. INTRODUCTION 

Narrow band surveys have been used to discover 
Lyq-emitting (LAE) galaxies at high r edshift (e.g., 
Hu fc M cMah onl lilM ICowie fc Hul fl998t ISteidel et al.1 
20001 iRhoads et all 120031 ) and to study their proper- 
ties. The galaxies' strong emission lines reveal a set 
of young, potential l y dus t-free galaxies theorized by 
iPartridge fc Peebles! (|1967t) . LAEs allow us to study 
galaxy formation beginning with one of the smallest 
building blocks found so far. At z ~ 3.1, typical LAEs 
are low mass galaxies with MstMar — 10 9 Mq and little 
dust extinction, Ay < 0.2 (|Gawisei1l2009t iNilsson et al.l 
l2007t lAcquaviva et al.ll2011al). These objects have been 
observed at z > 3 (iVenemans et al.ll2005t iGawiser et al.l 
l2006tlGronwall et al.ll2007t INilsson et al.ll2007t ) and as far 
out as z ~ 7 (jive et al.ll2006ir ~ 

Previous papers have studied the evolution of LAEs 
and other, gener ally more massiv e, high-z galax y pop- 
ulations (e.s. iGilh et al.1 I2007L iQuadri et al.l 120071 : 
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iSalvadori et al.l [20101 : lYaiima et al.ll2011l ) by using clus- 
tering properties as a technique to make evolutionary 
claims between redshifts. The connection between z = 
3.1 LAEs and prese nt-day galaxies was determined by 
IGawiser et all (|2007t ) (hereafter Ga07) by measuring the 
clustering properties of t h e z = 3.1 LAEs from the 
sample of iGronwall et all (120071 ) Ga 07 used the for- 
malism devised by IMo fc White! ()1996[ ) to compute the 
median dark matter mass of their host halos to be 
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IQuchi et al. (2010). A weakness of the analytical con- 
ditional mass function is its inability to distinguish in- 
dividual substructures within a halo or to predict the 
masses of these substructures. N-body simulations fill 
this gap by producing halo merger trees that allow us to 
determine properties for the descendant halos and their 
substructures. 

Spectral energy distribution fitting for z = 3.1 LAEs 
reveals a large range in stellar ages (G a07; lLai et all 
120081 INilsson et al.l [20071 : lOno et alJl2010T ). Ga07 deter- 
mined the young stellar component to have an age of 
20^ 3 q Myr (Ga07) using a two-stellar population model. 
This could have two simple interpretations: 1) LAEs oc- 
cur as a galaxy's first burst of star formation, which lasts 
~ 40 Myr, or 2) LAEs are a recurring phase, where each 
burst of star formation lasts ~ 40 Myr (Ga07). In both 
cases, stellar evolution produces dust, which ultimately 
quenches Lya emission. The old stellar population in 
the two stellar population model is not well constrained, 
with an age of up to 2 Gyr. A single stellar population 
model, which observations can not dism iss, has an age of 
50 — 1 Mvr (TLai et al.l 120081: s ee also lAcquaviva et al.l 
l2011a|) . lAcquaviva et al.l (l2011bD determined stell ar ages 
of ~ 1 Gyr for z ~ 3.1 LAEs. INilsson et al.l (pOOl found 
z ~ 3.1 LAE stellar ages of 0.85±g;^| Gyr. 

These clustering and spectral energy distribution re- 
sults allow us to create mo ck LAE catalogs within the 
Millennium-II SimulatiorE] (|Bovlan-Kolchin et al.l I2009L 
hereafter MS-II) to study the dark matter mass evolution 
to the present day within a large cosm o logical simulation. 
A similar investigation bv iCoil et al.l (|2008a! ) used star- 
forming galaxies at z ~ 2 as a starting point and studied 
their dark matter halo evolution until the present-day 
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within the Millennium simulation. Section 12.11 describes 
the MS-II. Section [3] presents specifics of the mock cat- 
alogs, while sections [4] and [5] report on the clustering 
analysis and descendants of the catalogs respectively. All 
distances reported are comoving, Ho = 100/ikm/s/Mpc 
and h = 0.73 is assumed throughout. 

2. SIMULATION AND ANALYSIS 

2.1. Millennium-II Simulation 

Our study uses the results of the MS-II run by the 
Virgo Consortium. The MS-II contains 2160 3 parti- 
cles in a cube of 100 h~ 1 Mpc on a side. The particle 
mass is 6.9 x 10 6 /i _1 Mq with a minimum halo mass of 
1.38 x 10 8 /i" 1 M Q . MS-II gives us the ability to resolve 
125 times less massive sub-halos than those observed 
within the Millennium Simulation (jSpringel et all 120051 : 
iLemson fc Virgo Consortium! 120061 ) . LAEs at z = 3.1 
appear to be hosted within low mass halos, making the 
improved mass resolution necessary. Although the Mil- 
lennium Simulation provides 125 times larger volume, 
MS-II offers robust statistics for spatial clustering of ha- 
los in the mass range of interest. Both simulations use a 
ACDM cosmology with values: 



Q tot = i.o, vt m = 0.25, n b = 0.045, n A 

h = 0.73, er 8 = 0.9, n s = 1. 
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All values except n s and as are within l a of the val- 
ues r eported in the 7- year WMAP results (jLarson et al.l 
MM- 

MS-II also offers improved temporal resolution, with 
dark matter halos selected in 67 timesteps using a 
friends-of-friends algorith m dDavis et al.l I1985T) with a 
linking length of b = 0.2 (jBovlan-Kolchin et al.l [20091) . 
Each FOF group was analyzed for sub-ha los using the 
SUBFIND algorithm (jSpringel et al.l 120011 ). which iden- 
tifies gravitationally bound sub-halos within the FOF 
group. We defined a central sub-halo to be the most 
massive substructure within an FOF group, and other 
sub-halos within the group were classified as satellites. 

During our age selection (described in £|3.2[) we found 
that some of the youngest dark matter halos exhibited 
rapid mass growth in the range of a factor of 10-1000 
during one timestep (approx. 200 Myr). We believe this 
mass growth is unrealistic and is caused by misidentifi- 
cation of ownership of dark matter particles by the SUB- 
FIND algorithm between neighboring sub-halos. The er- 
rant merger trees appear to normalize by merging with 
the mass theft victim after a couple of timesteps. This 
allows us to trust the results from our descendants that 
have had a few timesteps to normalize. Our selection 
method, that removes fast growing merger trees, com- 
putes the dark matter mass ratios between the z = 3.1 
dark matter halo and its most massive predecessor and 
between the most massive z = 3.1 predecessor and its 
most massive predecessor. If either of these ratios are 
greater than 10, then we remove this merger tree from 

1 Data created from the MS-II can be accessed from 

the Max-Planck Institute and Durham University servers 

at http://www.mpa-garching.mpg.de/galform/millennium-II us- 
ing a Structured Query language [SQL") query. For a 
merger tree explanation see http://www.g-vo.org/Millennium- 
II /Help?page=mergertrees . 



our catalogs. With this filtered set of dark matter halos 
we produce our models. 

2.2. Sub-Halo Abundance Matching 

We use Sub-Halo Abundance Matching (SHAM) to de- 
termine stellar masses for the descendants of our models. 
The main principle of the SHAM algorithm is to assign 
luminosities or stellar masses from a luminosity or stellar 
mass function to dark matter masses from N- body dark 
matt er simulations monotonically (e.g., iConroy et al.l 
120061 ) . Dark matter sub-halos are assigned stellar masses 
by matching the number densities of halos such that 



n h (> M DM ) = U g {> M S teUar(M DM )). 



(1) 



The dark matter masses used in eqn. [T] are sub-halo 
masses. If the sub-halo is a satellite substructure then 
we modify the dark matter mass to be the infall mass. 
Infall mass is defined as the larger of the mass of a cen- 
tral sub-halo before it falls into a larger halo to become 
a satellit e or the satellite ma ss in the following timestep 
CSee e.g.. lConrov et al.ll2006l ). We use the infall mass be- 
cause it gives a better representation of the current stel- 
lar component. Unlike the dark matter particles that are 
disrupted during infall, the stellar component is deeper 
in the gravitational potential well and therefore is less 
likely to be disrupted during infall. 

Our algorithm solves eqn. [T] by using a Newton's 
method approach where the function and its derivative 
are 
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(3) 



In the above we have used a Schecter function with 
parameters 0» , M* and a and T is the incomplete gamma 
function. This method allows us to assign stellar masses 
by solving eqn. [2] for Ms te ;; ar . The large number of sub- 
halos in the simulation forces us to use a sub-selection 
technique to determine the stellar mass to dark matter 
sub-halo mass function. Afterwards we interpolate this 
function to determine the stellar masses of the z — 
descendants from their sub-halo dark matter masses. 

3. LAE MOCK CATALOGS 

We created ten LAE mock catalogs from the MS-II 
using sub-halos at z = 3.06 (snapnuii031) based on mass 
and age selection. The catalog names and properties are 
listed in Tables Q] and [2j In the following sections we 
will discuss the selection criteria and motivation for the 
different catalogs. 

3.1. Mass-based Catalogs 
2 Snapnum is the snapshot (timestep) number from the MS-II. 
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All mass-based catalogs were chosen to reproduce the 
z = 3.1 LAE number density ( 1.5 x 10~ 3 Mpc~ 3 ) deter- 
mined bv lGronwall et all (|2007f ). thus generating models 
with 3856 sub-halos in the MS-II volume. 

1 . The mass limit criterion was selected to have a min- 
imum mass of log 10 M Min = 10.6 M to match that 
inferred by Ga07 from the observed clustering of 
z — 3.1 LAEs. This approa ch is similar to the one 
used by iCoil et all (|200 8a) ; their mass limit was 
chosen to reproduce the clustering of z ~ 2 star- 
forming galaxies. We randomly selected 5.24% of 
these sub-halos to match the above number density. 
The Mass limit catalog was designed to reproduce 
the observed correlation length of z = 3.1 LAEs. 

2. The Median catalog was selected using the median 
mass reported by Ga07, log 10 M Mcd = 10.9 M Q , as 
the catalog's median mass value. We expanded 
evenly around the median mass to obtain the ob- 
served LAE number density, choosing halos with a 
range of 7.56 x 10 10 M -8.36 x 1O 1O M . As ex- 
pected, this catalog and the Mass Limit catalog 
have similar median masses; it should also repro- 
duce the observed correlation length. 

3. The —a and +a catalogs were selected and named 
based on the ±ltr uncertainty reported by Ga07 
in the observed LAE median mass. The median 
masses of the catalogs are log 10 M Mcd = 10.0 M Q 
and 11.4 Mq, respectively. We expand evenly 
around their respective median masses to obtain 
the observed LAE number density. The mass 
ranges for the —a and +<r catalogs are 9.93 x 
10 9 M Q -1.01 x 10 10 M and 2.15 x 10 11 M Q -3.04 x 
10 11 M , respectively. These catalogs were created 
to study the uncertainty in the descendant prop- 
erties propagated from the observed clustering un- 
certainties and are not expected to reproduce the 
best- fit observed correlation length at z = 3.1. 

3.2. Age-based Catalogs 

Three age definitions were chosen to study the depen- 
dence of descendants' properties on age selection. The 
age definitions use merger trees rooted in central sub- 
halos within FOF halos with M FOF > 3.98 x 10 10 Mq, the 
mass limit criteria. Figure [T] shows an example merger 
tree showing the formation, assembly and merger ages 
assigned. Ages are defined as the difference of lookback 
time chosen and the lookback time at z = 3.1. The three 
age definitions are as follows: 

1. The Formation age (|Gao et alJl2004L I2005D quanti- 
fies the timescale for mass growth of the most mas- 
sive dark matter structure. It is assigned by finding 
the most recent timestep where the sub-halo mass 
of the most massive progenitor is less than half of 
the maximum mass in the entire merger tree. If 
this occurs between two timesteps, we linearly in- 
terpolate between these two to estimate the time 
when half the mass was accreted. 

2. The Assembly age (|Navarro "eTaIIIT997l) measures 
when half the maximum mass of a galaxy is present 
in collapsed sub-halos even if they have not yet 



merged. This is assigned by finding the most re- 
cent timestep where the sum of progenitor sub- halo 
masses is less than half of the maximum mass in 
the entire merger tree. If this occurs between two 
timesteps, we linearly interpolate between these 
two to estimate the time when half the mass has 
assembled. Because of these definitions, the As- 
sembly ages are always equal to or greater than 
the Formation ages. Figure [2] a) shows the relation 
between Formation and Assembly ages. 

3. Merger age searches for the most recent major 
merger. Our definition was designed to to use a 
main sub-halo to follow the m ajor merger which i s 
similar to the methods used bv lGenel et all ((2008). 
We define a major merger to occur when two cen- 
tral sub-halos have a mass ratio of 3:1 or less in a 
timestep and the most massive halo within that 
timestep is involved. The major merger begins 
when one of the descendants in the next timestep 
is a satellite of the other descendant or both cen- 
trals merge into a single sub-halo. In Figure [1] we 
have two central sub-halos, around 1.2 Gyr, where 
the lower mass central sub-halo descends into a 
satellite belonging to the descendant of the more 
massive central sub-halo. We track the descen- 
dants until the two sub-halos merge, possibly trig- 
gering star formation. We average the timestep 
of this sub-halo merger and the previous timestep 
to assign the merger time. The infall timescale 
from the beginning of the 3:1 merger of the two 
centrals to the final merger of their descendants 
is comp arable to the dynam ical friction timescale 
used by iGenel et all (120081) based on work from 
iBovlan-Kolchin et all (|2008D . Figure b) fc c) 
show little to no correlation between Merger age 
and Assembly or Formation ages. 

The age distributions are presented in Figure [3] a). The 
LAE mock catalogs are then created from the ages calcu- 
lated for all halos with the mass limit criteria by selecting 
the 5.24% youngest and 5.24% closest to median aged 
sub-halos. See Figure [3] b), c) & d) for age distributions 
for all halos and the age distributions of the selected cat- 
alogs. Table [2] lists the age properties of the different 
catalogs with the bold values being age statistics used to 
select that particular catalog. 

All the median age-selected catalogs have median ages 
of ~ 1 Gyr. For star formation triggered by accretion 
of sub-halos in minor mergers, Formation and Assem- 
bly ages should roughly track the age of stars in a 
galaxy i. e. the population domina ting its stellar mass 
(see e.g., iDe Lucia fc Blaizotll2007l ). Since SED fits al- 
low ~ 1 Gyr old populations t o comprise the majority 
of stellar mass in LAEs (Ga07: lAcauaviva et alj|2011bt 
iNilsson et all 12007). the median Formation and Assem- 
bly models are feasible models. However, because the 
merger age should track a young stellar population born 
in a starburst triggered by the major merger, we can- 
not reconcile a ~ 1 Gyr median merger age with SED 
results; this rules out the Median Merger model at high 
confidence. The young Formation and Assembly catalogs 
are viable models which have median ages of 240 Myr, 
consistent with ages found for single-population models. 
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The young Merger catalog has a median age of 91 Myr 
and is consistent with the observed age of the young stel- 
lar component. 

4. CLUSTERING ANALYSIS OF MOCK CATALOGS 

We used the naive estimator £jv (e.g. JLandv fc Szalavl 
I1993T ) to calculate the two point auto-correlation function 
(2PCF) for the models. 



6v(0 



DD(r) 
RR{r) 



DD(r) = 
RR(r) = 



2 X Number of data-data pairs within a radial bin 

noino - 1) 
Volume contained within a radial bin 

Volume of simulation 

Data-data pairs are unique pairs of sub-halos, from a 
catalog, that arc separated and binned by radius r. The 
simple geometry of the simulation, in conjunction with 
the periodic boundary conditions, allow us to use a geo- 
metrical formula for RR, which eliminates uncertainty in 
the estimator caused by binning random-random pairs. 

We applied a correction to our data abscissas to match 
the effective center of the radial bin using 



bin ~ 
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The correction was applied by our fitting algorithm to 
shift the data point to the radius which corresponds to 
the average value of the bin. We fit our naive estimator 
with the power law given by 



(r/ro) 7 - ljq 



1 + wn 



RR{r) 



ro 



dr. 



(6) 
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where is the integral constraint found from the power- 
law term, (7^) 7 for the parameters during fitting. We 
minimized \ 2 to determine the best fit for the parameters 
ro and 7 using 



1 - RR{r) 



(8) 



l + LJa ) \{n D {n D -l)/2)RR(r) 

for th e variance of the naive estimator (|Landv fe Szalavl 
[1991 . 

In the following discussion we will fix 7 = 1.8 to com- 
pare our results to the observed correlation length for 
z = 3.1 LAEs, r = 3.6^ Mpc (Ga07). Table Q] lists 
the values found for ro fixing 7 = 1.8 and also fitting 
both ro and 7 as free parameters, while Figures 2] and [5] 
show the best fit 2PCF for the catalogs. 

4.1. Clustering of Mass-Based Catalogs 

Figure @] and Table [TJ show the mass-based catalogs' 
2PCFs, correlation lengths, 7 and x 2 values. We con- 
firm that correlation lengths increase with median mass 



in our mass-based catalogs for both fixed 7 = 1.8 and 
when 7 is allowed to be a free parameter. The Median 
and Mass limit catalogs have similar median masses and 
were expected to have the same correlation lengths when 
7 = 1.8; the correlation lengths are consistent with one 
another. These two models also have correlation lengths 
that are consistent with the observed LAE correlation 
length of r = 3.6T. X ; Mpc, making them good represen- 
tations of LAEs at z = 3.1. The +<r and— a models have 
correlation lengths close to the ±1<7 uncertainty from ob- 
served values, as expected. Fixing 7 = 1.8 allows us to 
compare our results to previous works, but the best fit 
models are not consistent with this value. We find an 
average val ue of 7 = 1.33, which i s con sistent with the 
results from lBovlan -Kolchi n et all ([2009) at z = 2.07 us- 
ing our fitting range. 

4.2. Clustering of Age-Based Catalogs 

The similarities of the mass properties of the age-based 
catalogs (see Table Q} imply that any difference in the 
clustering properties between age-based catalogs are due 
to the age definition and selection. The Formation Me- 
dian and Assembly Median models have similar correla- 
tion lengths of ro ~ 3.6 Mpc and are consistent with the 
observed LAE correlation length. The Assembly Young 
and Formation Young have higher correlation lengths of 
ro ~ 4.4 Mpc and are barely consistent at the 68% level 
with the z ~ 3.1 LAE correlation length. We cannot rule 
out the Young Merger model at 95% confidence. The Me- 
dian Merger model has a correlation length large enough 
to be ruled out at 93% confidence. This large correla- 
tion length and the high merger ages make the Median 
Merger model an unrealistic z ~ 3.1 LAE model. No 
other models are ruled out based on clustering. The ob- 
served trend for the Formation and Assembly ages that 
young aged models have a larger clustering length com- 
pared to the m edian aged model s does not disagree with 
the findings of iGao et al.l ()2005l ) that the youngest 20% 
of halos at z — cluster less strongly than the oldest 
20%. Once we use a similar selection we also find that 
the oldest halos have higher clustering compared to the 
youngest. 

5. DESCENDANTS OF MOCK LAES 

We used the MS-II merger trees to find the z = 2.07 
(snapnum 36), z — 0.99 (snapnum 45), and z = (snap- 
num 67) descendants of the mock LAE catalogs. We 
report the mass distributions of the descendants based 
on the mass of the host FOF group, M FOF . We prefer 
M FO f because it determines the evolution of dark matter 
halos and is used in the Press-Schechter formalism, while 
the individual sub-halo mass traces the evolution of in- 
dividual galaxies. We classify the most massive sub-halo 
within each FOF group as a central and other sub-halos 
as satellites. All other smaller sub-halos within the FOF 
group are satellites. The following sub-sections will dis- 
cuss the descendants of our mock catalogs as reported in 
Tables Hi 

5.1. Mass Limit and Median Mass Model Descendants 

Figure [5] shows an example histograms demonstrating 
the evolution of the Mass Limit and Median models from 
z — 3.1 to z — using FOF halo masses. Both models 
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have the same median masses, logjn M me( ; ~ 10.9 M©, 
at z = 3.1, as also seen in Table [1] All other models 
evolve in a similar manner. Figure [Jj (bottom panel) 
summarizes this same evolution for all of the mass se- 
lected models. The models' descendants also have similar 
satellite mass distributions (dashed histograms in figure 
[5]) as they evolve with time. The satellite median masses 
grow from M FOF . Mcd = 10 12 M© at z = 2 to 10 13 - 7 M© at 
z = 0. The central descendants show evolution towards 
higher masses, as expected for bottom-up halo growth 
(|Davis et al.l 119851 ). though the 10th percentile in mass 
stays roughly constant at its z ~ 3.1 value. We find that 
the catalogs' central population, which comprise 55-60% 
of the descendants, has mass growth from M FOF; Mcd — 
10 10 - 9 M© at z = 3.1 to Milky Way-sized dark matter 
halos, ~ 10 n - 8 M©, at z = 0. The median FOF mass 
for all descendants grows from M FO f ; Mcd 

= 10 10 - 9 M© at 
z = 3.1 to MpQF; Mod = 10 12 - 6 M© at z = 0. 

For comparison, Figure [7] top panel shows these de- 
scendants in terms of the individual sub-halo mass in- 
stead of the mass of the FOF group that the sub-halo 
resides in. The evolution of mass in this sense is similar 
to the FOF mass evolution; sub-halos tend to grow in 
mass towards z = 0, but mass loss is seen in the satellite 
populations. Central sub-halo masses are similar to the 
FOF masses; therefore, we expect a similar median mass 
at z = 0. We find the central sub- halo median mass to 
be M = 10 n - 8 M© at z = 0. Comparing the top and bot- 
tom panels of figure [7] we see that descendants that are 
satellites reside in massive halos, but the satellite sub- 
halos are less massive, with a median mass of 10 108 M© 
at z = 0. 

5.2. Age-Based Model Descendants 

The similarities in the six age-based models' mass dis- 
tributions at z = 3.1 suggests that we should expect 
similar descendant distributions. Figure [5] top and bot- 
tom panel confirm the similarities among these catalogs' 
descendants. The mass evolution represented in figure |5] 
bottom panel is the mass of the FOF halo where the sub- 
halo resides while the top panel uses sub-structure mass. 
The evolution of the mass distributions are all similar and 
therefore independent of age definition or selection. The 
descendant central median masses are M FOF = 10 112 M© 
at z = 2.1, 10 n - 5 M© at z = 1, and 10 12 M© at z = 
with a small scatter of less than a factor two. We find 
that the fraction of central sub-halos for all age-based 
models is independent of age definition or selection (Ta- 
bles |3][S]) . The satellite descendants at z = comprise 
42 — 44% of the population and have a median FOF mass 
of 10 137 M©. For all the models the full distribution of 
descendants have a median FOF mass within a factor of 
two of 10 12 7 M© at z = 0. These FOF masses are in 
agreement with those found for the Median and Mass 
limit catalogs described in section 15.11 The sub-halo 
masses of the age-based models are also similar to the 
values obtained for the Median and Mass limit models. 

5.3. +a and —a Model Descendants 

The +<7 and —a models were created to quantify the 
uncertainties within the descendant distributions caused 
by the uncertainties in the observed clustering analysis. 
Figure [Jj shows the +a~ and —a with the other mass- 



selected models for comparison using the FOF mass. The 
central descendant sub-halos have median masses that 
are a factor of four higher for the +a model and a factor 
of ten lower for the —a model compared with the Median 
mass catalog's descendants at z = 0. The satellite dis- 
tributions at z = have median masses that are a factor 
of three larger for the +a and a factor of five smaller for 
the —a than the Median model. When we consider the 
full descendant distribution the median mass of the +a 
model is a factor of three larger and the —a model is 
six times lower than the Median model. The fraction of 
central descendant sub-halos is similar between the two 
models. 

6. DETERMINING STELLAR PROPERTIES OF 
DESCENDANTS 

Using the SHAM algorithm discussed in section 12.21 
we determine the stellar masses for all sub-halos at z = 
0. We use the Schecter funct ion parameters fou nd by 
the Sloan Digital Sky Survey (jPanter et all 120071 ). The 
parameters are 

0* = 2.2 ± 0.5 stat ± l sys x 10- 3 Mpc- 3 

M* = 1.005 ± 0.004 stat ± 0.200 sys x 10 11 M© 

a = -1.222±0.002 stot ±0.1 SJ/s . 

Figure O shows the relationship between infall and stellar 
masses determined from the SHAM algorithm. 

After application of the SHAM algorithm to the entire 
z = MSII sub-halo catalog covering dark matter masses 
ranging from 10 8 28 M© to 10 14 94 M©, we obtained stel- 
lar masses ranging from 0.14 M© to 10 1180 M©. The very 
low stellar masses are attributable to the excess of dark 
matter sub-structures at low mass compared to the num- 
ber of galaxies from the stellar mass function, which re- 
mains an unsolved problem in galaxy formation. The 

— la model is the only model with a significant number 
of halos in this range. Figure [TU] shows the median stellar 
masses for all the models. The error bars denote the 10th 
and 90th percentiles. All the models except the +la and 

— la have similar stellar mass distributions. The +la 
and —la models have corresponding shifts in their me- 
dian stellar masses due to their selection. 

7. DISCUSSION AND CONCLUSIONS 

We find that the Median, Mass Limit, Young Forma- 
tion, Median Formation, Young Assembly and Median 
Assembly models have correlation lengths within the 68% 
confidence interval for LAEs at z — 3.1, all though the 
Young Formation and Young Assembly models have cor- 
relation lengths near the +la limit. As designed, the +a 
and —a catalogs have z — 3.1 correlation lengths near the 
observed ±ltr correlation length uncertainty. The Young 
Merger and Median Merger models have large correlation 
lengths and are not consistent with the z — 3.1 clustering 
measurement. We eliminated the Merger Median model 
due to the age of ~ 1 Gyr at z = 3.1, in contrast with 
starburst ages of 20-100 Myr determined from SED fit- 
ting. However, SED modeling has not yet constrained 
the age of LAEs sufficiently to rule out the other age- 
based catalogs. 

We studied the connection between LAEs at z — 3.1 
and galaxy populations at z — 2.1 for which previous 
studies have determined typical dark matter halo masses. 
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iGuaita et all ()20100 found LAEs at z = 2.1 to have a me- 
dian dark matter halo mass of log(M/M Q ) = 11.5~^ 5 . 
Their study used bias evolution from the conditional 
mass function to show that z = 2.1 LAEs are possible 
descendants of z — 3.1 LAEs observed by Ga07. All of 
our models, except the +<r and —a models, have median 
dark matter halo masses of the full descendant distribu- 
tions around 10 n - 3 M Q at z = 2.1. Hence, z = 2.1 LAEs 
could be direct descendants of z — 3.1 LAEs, rather than 
a new set of halos undergoing their first phase of star for - 
mati on. However, S EP m odeling bv IGuaita et al.1 (|2011f ) 
and iNilsson et al.l (|2011[ ) have found typical starburst 
ages for z = 2.1 LAEs of 12I^ 49 Myr and 80^0 M Y r - 
Given these starburst ages, the z = 2.1 LAEs would have 
to be experiencing a subsequent burst of star formation 
to be descendants of z = 3.1 LAEs. Another study by 
( Adcl berger et al.ll2005l see their Fig. 2) determined the 
clustering of BX galaxies at z ~ 2. The clustering re- 
sult for the least luminous (K s > 21.5) subset implies a 
median dark matter mass of log 10 M Mod /M Q = ll.Oj^g. 
This result is consistent with our z = 2.1 descendant me- 
dian dark matter mass, except for the +er and —a mod- 
els. We find that ~ 75% of the descendants of z = 3.1 
LAEs are centrals at z = 2.1; it is unclear whether ~ 2 
LAEs and BX galaxies represent a similar mix of centrals 
and satellites. 

We can make a similar comparison from z — 2.1 to 
z = 0.9. The Deep Extrag alactic Evolutionary Probe 2 
(DEEP2. lDavis et aT1l2003f) measured the clustering bias 
of color selected galaxi es at z ~ . 9. For our chosen cos- 
mology, their results ([Coil et al.l l2008bf) imply median 
dark matter halo masses of M Mcd = io 12 - 9±0 - 1 M Q and 
M Med = 10 12 0±0 - 1 M for red and blue galaxies, respec- 
tively. The Mass limit, Median, and age-based models 
have z = 1 descendant median dark matter halo masses 
of M ~ 10 119 M Q , consistent with the blue galaxies from 
the DEEP2 survey. The DEEP2 red galaxy subset are 
not consistent with being descendants of our LAE cat- 
alogs due to their larger dark matter mass. LAEs at 
z = 3.1 are therefore possible progenitors of blue (late- 
type) galaxies residing in dark matter halos with mass 
1O 12 M at z ~ 1. We find that - 65% of the descen- 
dants of z = 3.1 LAEs are centrals at z = 1; it is unclear 
whether the blue DEEP 2 galaxies represent a similar mix 
of centrals and satellites. 

The models that best represent LAEs at z = 3.1 pro- 
duce z = descendants with a median mass around 
10 12 7 M Q . We find that ~ 55% of the descendants at 
z = are central sub-halos, with little variation on the 
fraction of the central descendants based on model selec- 
tion. When we study only centrals, the age-based mod- 
els have median dark matter halo masses at z — of 
10 12 M while the Median Mass and Mass limit catalogs 
have slightly lower values of 10 118 M©. We see no sig- 
nificant dependence on age definition or selection of the 
descendants' mass distributions. These results show that 
the Mass Limit, Median Mass, and age-based catalogs' 
descendants have a majority of central sub-halos which 
reside in L* type dark matter halos at z = 0. How- 
ever, the +er and —a models have descendant masses 
that are a factor of three larger and seven times smaller 
than the median mass of the Median model. The factor 
of five spread in median mass of the +a and —a cata- 



logs' present-day descendants shows that a more precise 
measurement of the z ~ 3 LAE bias is needed to place a 
stronger constraint on the descendant properties. 
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Halo Mass 



Figure 1. Representation of a merger tree which shows the ages found for one particular z = 3.1 central sub-halo. Sub-halos are 
represented by circles, where filled circles are centrals and open circles are satellites. The black, blue and red horizontal dashed lines mark 
the Formation, Assembly and Merger ages respectively. The black vertical dashed line marks half the maximum mass of sub-halos in the 
merger tree. The Formation age is assigned when the most massive sub-halo crosses this vertical line. The blue solid curve is the sum of 
the sub-halo masses at each timestep and the Assembly age is marked by the intersection of this curve with the black vertical line. The 
red solid curve marks the left boundary of the region where a 3:1 ratio between the most massive object and any other sub-halo can occur. 
The sub-halos involved in the beginning of the < 3 : 1 major merger are marked with red circles. The Merger age is assigned when the 
< 3 : 1 ratio satellite disappears due to having fully merged with the central. 




Figure 2. Figures show the relationships between the different age definitions. Red curves are the median with 25th and 75th percentiles 
shown as red error bars. Panel a: Shows the relation between the Assembly and Formation ages. By definition the Formation age is equal 
to or greater than the Assembly age. Panel b: Shows the relation between the Assembly and Merger ages. There is no correlation between 
these two age definitions. Panel c: Shows the relation between Formation and Merger a ges. There is little to no correlation between these 
two age definitions. The merger ages are discrctized due to averaging of timesteps (See H3.2I I. 
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Figure 3. Light green, peach and royal blue curves are the age distributions for central sub-halos at z = 3.1 with M f of > 3.98 X 10 10 
Mq using the Formation, Assembly and Merger definitions. Panel a: Shows the three age distributions. Panels b,c,d show the Formation, 
Assembly and Merger age distributions (same colors as panel a) along with those of the selected young and median aged mock halo models 
in different colors. 




1 101 10 

Radius (Mpc) 

Figure 4. The two-point autocorrelation functions for mass-based mock models. See Table [l] for correlation length and gamma values. 
Error bars are a function of the model being fitted. Panel a): We fitted a power-law with 7 = 1.8 to the data to determine the correlation 
length for each model. As expected, the models show a trend to larger correlation length with higher median mass. Panel b): We fit a 
general power-law to determine the best fit correlation length and 7. We find 7 ~ 1.4 for these mock LAE models. 
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Figure 5. Two-point autocorrelation functions for age-based catalogs. See Table [T] for correlation length and gamma values. Panels a, 
b, and c): We fit a power-law with 7 = 1.8 to the data to determine the correlation length for each model. The correlation length of the 
formation and assembly age models separate based on age selection. The merger age definition shows a smaller separation based on age 
selection. Panels d, e. and f): We fit a general power-law to determine the best fit correlation length and 7. We find 7 ~ 1.4 for these LAE 
catalogs. We find an age-clustering relation when we use the formation and assembly ages where the youngest dark matter halos are more 
clustered than the median aged dark matter halos. 
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Figure 6. Mass histogram, using FOF mass, for the Median (dark red) and Mass limit (green) mock catalogs at z = 3.1 (upper-left) and 
descendants at z = 2.1 (upper-right), 2 = 1 (lower-left) and z = (lower-right). Solid and dashed curves represent the central and satellite 
mass distribution breakdown. Satellite sub- halos have a larger mass than central sub-halos because we plot MpQp, the mass of the FOF 
group to which the sub-halo belongs. Vertical dashed lines delineate the median mass of each distribution. 
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Figure 7. Figure shows the sub-halo (top panel) and FOF (bottom panel) mass evolution of mass selected models through redshift. The 
Median (dark red), Mass limit (green), +lcr (orange) and —la (black) models are evenly spaced around their redshift for easier viewing. 
Solid and open points represent the distributions for centrals and satellites, respectively. Points are the median sub-halo masses while the 
error bars show the 10th and 90th percentiles for the distributions. The sub-halo and FOF masses for centrals tend to increase towards 
lower redshift as sub-halos merge and accrete. Some of the z = 3.1 central halos merge with more massive halos becoming satellites. These 
satellites reside in massive FOF halos, but have a smaller individual sub-halo mass. 
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Figure 8. Figure shows the sub-halo (top panel) and FOF (bottom panel) mass evolution for our age selected models. All the age selected 
models have similar z = 3.1 masses. Median sub-halo mass values tend to increase towards low redshift as the dark matter halos merge 
and accrete. There is a small difference in the median masses between the young formation/assembly and median formation/assembly due 
to age selection. This small mass difference propagates to z = in the median masses of the centrals. Otherwise these models have nearly 
the same distribution at each redshift. Some of the z = 3.1 central halos merge into more massive FOF halos producing satellites at lower 
redshift. These satellites have smaller sub-halo masses due to merging into a more massive FOF halo. 
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Figure 9. Relation between infall and stellar mass found from the SHAM algorithm after application to all z = sub-halos. We use this 
relation to determine the stellar ma sses from the z = descendant sub-halos. The stellar mass function used in the SHAM algorithm had a 
fitting range of 10 8 - 5 - 10 1185 M l|Panter et al.ll2007ft and stellar masses below and above these values are extrapolations from the SDSS 
fit. 
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Figure 10. Median stellar masses of the descendants at z = 0. The error bars denote the 10th and 90th percentile stellar masses for each 
of the models' distributions. All models except the + 1ct and — la models have similar stellar mass distributions. The low stellar masses 
found in the — l<j model are caused by the extrapolation of the stellar mass function towards low mass (See JHJ. 
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Table 1 

Mass Properties and correlation lengths for z = 3.1 LAE models 



Name a 


25th % Mass b 


Median Mass b 


75th % Mass b 


ro c 


7 C 


X 2 /D.O.F. 


r d 


7 d 


X 2 /D.O.F. 




log lo (M ) 


log 10 (M Q ) 


log 10 (M Q ) 
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Table 2 

Age Properties of z = 3.1 LAE models 



Name" 


Form. 25th 
Gyr 


Form. 50th 
Gyr 


Form. 75th 
Gyr 


Asscm. 25th 
Gyr 


Asscm. 50th 
Gyr 


Asscm. 75th 
Gyr 


Mcrg. 25th 
Gyr 


Mcrg. 50th 
Gyr 


Merg. 75th 
Gyr 


+lv 


0.523 


0.696 


0.857 


0.644 


0.864 


1.033 


0.583 


0.978 


1.454 


Mass Limit 


0.569 


0.748 


0.921 


0.724 


0.946 


1.110 


0.583 


1.090 


1.527 


Median Mass 


0.580 


0.766 


0.942 


0.754 


0.964 


1.128 


0.583 


1.090 


1.527 


-lcr 


0.675 


0.877 


1.066 


0.915 


1.080 


1.218 


0.725 


1.288 


1.802 


Formation Age 




















Young Formation 


0.156 


0.258 


0.301 


0.164 


0.289 


0.354 


0.267 


0.856 


1.454 


Median Formation 


0.738 


0.747 


0.755 


0.831 


0.905 


1.018 


0.583 


0.978 


1.454 


Assembly Age 




















Young Assembly 


0.158 


0.270 


0.324 


0.164 


0.289 


0.345 


0.583 


1.090 


1.593 


Median Assembly 


0.684 


0.785 


0.847 


0.936 


0.946 


0.956 


0.583 


0.978 


1.454 


Merger Age 




















Young Merger 


0.419 


0.553 


0.686 


0.693 


0.897 


1.069 


0.091 


0.091 


0.091 


Merger Median 


0.640 


0.832 


0.981 


0.749 


0.970 


1.120 


0.978 


1.090 


1.090 



Note. — Tabic reports the 25th, 50th and 75th percentiles for each model using the Formation (Form.). Assembly (Asscm.) and Merger (Mcrg.) 
age definitions (Sec scction |3.2| l. Bold entries mark the values in the age definition used to select the age-based models. 
a Model descriptions found in Tabic fTl 



Table 3 

z = 2.1 Descendant Mass Distributions 



Name a 


Distribution Type b 


Fraction 


25th % Mass d Median Mass d 75th % Mass d 








logmfMm) 


log-m (M m ) 


-1-1 IT 


Full 


100% 


11.54 


11.66 
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10.24 




Qnfpllitp 


98% 


10.76 


11.24 


11.95 


Yn nncf T5Vi t* m atinn 

X V J LI J 1>A J- VJl 11 Id L> I VJ 11 


Full 


100% 


11.01 


11.36 


11.88 




C&TltTCLl 


72% 


10.93 


11.18 


11.59 




Qnfpllitp 


98% 


11.50 


11.91 


12.51 


1\/T prlian T5Vi t* m atinn 

± V J_ v J L 1 1 d. 1 1 J_ Ul 111CL L> JtJ JJ 


Full 


100% 


10.94 


11.25 


11.74 




C&TltTCLl 


77% 


10.88 


11.11 


11.43 




tint f>1 lit p 


23% 


1 1 fit; 
11. DO 


Iz.UO 


1 O C^O 

LZ.OZ 


Young Assembly 


Full 


100% 


11.01 


11.36 


11.91 




Central 


72% 


10.92 


11.19 


11.58 




Satellite 


28% 


11.52 


11.94 


12.55 


Median Assembly 


Full 


100% 


10.97 


11.26 


11.71 




Central 


77% 


10.91 


11.12 


11.44 




Satellite 


23% 


11.57 


11.97 


12.47 


Young Merger 


Full 


100% 


10.96 


11.30 


11.79 




Central 


75% 


10.90 


11.12 


11.48 




Satellite 


25% 


11.64 


12.04 


12.56 


Merger Median 


Full 


100% 


10.94 


11.24 


11.76 




Central 


76% 


10.89 


11.10 


11.40 




Satellite 


24% 


11.65 


12.06 


12.58 


a Model descriptions found 


in Tabic ^ 










b Centrals are the most m; 


issivc sub-halos within their FOF ha] 


o. Satellites arc 


all other sub-halos. 


Full describes the 


properties of all descendant 


sub-halos. The median 




for the central e 


Lnd satellite distribi 




Figures [7] and [8] 












C Fraction of objects within 


full, central or satellite 










Mass reported is the mas; 


3 of the FOF halo where 


the sub-halo 


resides. MpQp. 
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Table 4 

z = 1 Descendant Mass Distributions 



Name a 


Distribution Type° 


Fraction 


25th % Mass d Median Mass d 75th % Mass d 






logm (Mm) 


logmfMm) 


log-m (M m ) 


-1-1 IT 


Full 


100% 


11.83 


12.21 


12.90 




Central 


63% 


11.74 


11.92 


12.20 




^ntpllifp 


37% 


12.70 


13.10 


13.54 


1V/I sec Tiimit 


Full 


100% 


11.30 


11.86 


12.68 




t.pTt irn 1 

W O / Lb I Urh 


62% 


11.12 


11.43 


11.88 




kj UjLCL lilLC 


38% 


12.23 


12.81 


13.42 


\ T o r 1 1 '."i i \ iV/Tqcc 


Full 


100% 


11.27 


11.68 


12.53 






61% 


11.18 


11.34 


11.63 






39% 
oy /o 


12.21 


12.74 


13.33 


-la 


Full 


100% 


10.29 


10.70 


11.93 




C&TltTQX 


61% 


10.21 


10.34 


10.61 




kj LlliCL lilLC 


39% 
oy /o 


11.43 


12.18 


13.04 


Yn nncf T5Vi i* m atinn 

X V J LI J 1>A J- 11 Id L> V VJ 11 


Full 


100% 


11.40 


12.00 


12.79 






60% 


11.17 


11.54 


12.05 




kj UjLCL lilLC 


40% 


12.27 


12.80 


13.44 


TVT t> fi 1 £i n T5Vi t* m atinn 

± V J_v J L 1 1 CX. L 1 J_ Ul J J JCl L> JU J J 


Full 


100% 


11.30 


11.84 


12.63 




CcTltvoX 


63% 


11.13 


11.45 


11.89 




nipl Hip 


37% 




lz. to 


lo.oo 


Young Assembly 


Full 


100% 


11.40 


12.00 


12.79 




Central 


60% 


11.18 


11.55 


12.06 




Satellite 


40% 


12.26 


12.79 


13.44 


Median Assembly 


Full 


100% 


11.32 


11.83 


12.61 




Central 


63% 


11.15 


11.46 


11.87 




Satellite 


37% 


12.19 


12.71 


13.38 


Young Merger 


Full 


100% 


11.34 


11.91 


12.67 




Central 


62% 


11.16 


11.49 


11.97 




Satellite 


38% 


12.19 


12.75 


13.39 


Merger Median 


Full 


100% 


11.27 


11.83 


12.61 




Central 


63% 


11.12 


11.42 


11.87 




Satellite 


37% 


12.23 


12.78 


13.45 


a Model descriptions found 


in Tabic ^ 










As defined in Table 3. The median mass values foi 


: the central 


and satellite distr: 




in Figurcs[7]and[8] 


C Fraction of objects withir 


i full, central or satellite 


subset . 








^ Mass reported is the mas 


s of the FOF halo where 


the sub-halo 


resides, MpQp. 
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Table 5 

z = Descendant Mass Distributions 



Name a 


Distribution Type° 


Fraction 


25th % Mass d Median Mass d 75th % Mass d 






logm (Mm) 

o JLU \ l_J / 


l0K in (M(T,) 


log-m (M m ) 


-1-1 IT 


Full 


100% 


12.26 


12.95 


13.83 




Central 


60% 


12.05 


12.38 


12.93 




^ntpllifp 


40% 


13.41 


13.88 


14.36 


1V/I sec Tiimit 


Full 


100% 


11.75 


12.65 


13.69 




t.pTt irn 1 

W O / Lb I Urh 


58% 


11.42 


11.93 


12.61 




kj UjLCL lilLC 


49% 


13.04 


13.70 


14.32 


\ T o r 1 1 '."i i \ iV/Tqcc 


Full 


100% 


11.68 


12.51 


13.57 






55% 


11.44 


11.77 


12.32 






io /o 


12.97 


13.60 


14.13 


-1(7 


Full 


100% 


10.63 


11.71 


13.26 




C&TltTQX 


57% 


10.42 


10.73 


11.50 




kj LlliCL lilLC 


io /o 


12.29 


13.24 


13.96 


Yn nncf T5Vi t* m atinn 

X V J LI J 1>A J- 11 Id L> x yj 11 


Full 


100% 


11.91 


12.84 


13.74 






57% 


11.53 


12.07 


12.86 




kj UjLCL lilLC 




13.08 


13.70 


14.23 


TVT t> fi 1 £i n T5Vi t* m atinn 

± V J_v J L 1 1 CX. L 1 J_ Ul J J JCl L> JU J J 


Full 


100% 


11.79 


12.64 


13.66 




CcntvoX 


56% 


11.45 


11.94 


12.55 




1 aipl Hip 


44% 


i o no 


1 Q KQ 
lo.Do 


1 A OO 

14. zz 


Young Assembly 


Full 


100% 


11.92 


12.85 


13.76 




Central 


57% 


11.56 


12.08 


12.86 




Satellite 


43% 


13.09 


13.72 


14.32 


Median Assembly 


Full 


100% 


11.79 


12.60 


13.58 




Central 


58% 


11.47 


11.92 


12.52 




Satellite 


42% 


13.00 


13.63 


14.20 


Young Merger 


Full 


100% 


11.83 


12.71 


13.70 




Central 


58% 


11.47 


11.99 


12.67 




Satellite 


42% 


13.03 


13.71 


14.20 


Merger Median 


Full 


100% 


11.75 


12.56 


13.73 




Central 


58% 


11.42 


11.88 


12.51 




Satellite 


42% 


12.96 


13.77 


14.37 


a Model descriptions found 


in Tabic ^ 










As defined in Table 3. The median mass values foi 


: the central 


and satellite distr: 




in Figures [7] and [8] 


C Fraction of objects withir 


i full, central or satellite 


subset . 








^ Mass reported is the mas 


s of the FOF halo where 


the sub-halo 


resides, MpQp. 







